Time-resolved characterization of the formation of a collisionless shock 
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We report on the temporally and spatially resolved detection of the precursory stages that lead 
to the formation of an unmagnetized, supercritical collision-less shock in a laser-driven laboratory 
experiment. The measured evolution of the electrostatic potential associated with the shock unveils 
the transition from a current free double layer into a symmetric shock structure, stabilized by ion 
reflection at the shock front. Supported by a matching Particle-In-Cell simulation and theoretical 
considerations, we suggest that this process is analogeous to ion reflection at supercritical collisionless 
shocks in supernova remnants. 



Collision-less shocks represent particularly intriguing 
phenomena in plasma physics, due to their implica- 
tions in a broad range of physical scenarios, extending 
from laboratory-based laser-plasma experiments to as- 
trophysics. In the latter case, particular attention has 
been devoted to shock waves generated during the prop- 
agation of a supernova remnant (SNR) blast shell into 
the interstellar medium (ISM), since they are thought 
to be the dominant source of galactic high energy cos- 
mic rays [l|-|5[. In this case, non-collisionality is guaran- 
teed by the low particle collision frequency in the ISM 
@ ; the dynamics of SNR shocks is expected to be dom- 
inated by electromagnetic fields, thus setting stringent 
limits on their speed and stability. However, despite 
the considerable number of recent observations of such 
structures, the intrinsic difficulty in directly probing the 
plasma conditions around the SNRs has left the debate 
upon their generation mechanism and dynamics still open 
in the scientific community. A possible solution to this 
impasse might be provided by studying small-scale repro- 
ductions of these phenomena in laser-based laboratory 
experiments. Within this framework, promising results 
have indeed already been obtained 0-Q, especially con- 
cerning the stationary stage of these structures. Nonethe- 
less, a full understanding of the dynamics of a collision- 
less shock would also require a detailed characterization 
of the transient phase in which it is formed, a regime that 
has hitherto eluded experimental detection. 

Employing a time-resolved proton imaging technique 
[10( 1 , we present here the first experimental observation of 
the precursory stages that lead to the generation of a su- 
percritical electrostatic collision-less shock at the bound- 
ary of a blast shell of laser-ablated plasma expanding 
into a dilute ambient medium. By following the tempo- 
ral evolution of the propagation speed and of the pro- 
file of the associated electrostatic potential, it has been 
possible to distinguish the intermediate steps that let 
an initially freely expanding structure, similar to a cur- 



rent free double layer (CFDL) [Tl], [T2|, evolve into a 
forward shock propagating at a supercritical speed (ob- 
served Mach number of the order of 4). A matching 
Particle-In-Cell (PIC) simulation allowed us to reproduce 
the transformation of the plasma contact boundary by 
collision-less electrostatic processes into a forward and 
reverse shock on spatio-temporal scales that are similar 
to those in the experiment. The late-time behavior of 
such a shock, which is revealed by the simulation, shows 
that the shock-reflected ions locally heat the plasma in 
the close vicinity of the shock (forcshock) and thus in- 
crease the effective ion-acoustic speed, if compared to 
the one of the unperturbed upstream plasma. This heat- 
ing is not related to beam-driven instabilities but to the 
electrostatic potential associated to the shock-reflected 
ion beam and the increased ion-acoustic speed directly 
ahead of the shock contributes to its stabilization. This 
effect may also help stabilizing supernova remnant (for- 
ward) shocks like that of RCW86 [4J by heating up the 
ISM plasma directly ahead of the shock. 

In the experiment, carried out at the VULCAN laser 
facility at the Rutherford Appleton Laboratory [13], a 
1 ns, 120 J laser pulse (Long pulse in Fig. [T]) was fo- 
cussed, down to a focal spot of 50 /im radius, onto a 50 
/im thick gold foil, resulting in a peak intensity of the or- 
der of Tl ~ 10 15 Wcm -2 . The gold foil was surrounded 
by nitrogen gas at a pressure of 10 _1 mbar. A thin layer 
of hydrocarbon and vapor contaminants (thickness of the 
order of tens of nm) , as resulting from deposition of mois- 
ture and impurities of the low-density gas [Tj, [l5( , will be 
present on the gold foil surfaces. The laser-foil interac- 
tion generates a warm and dense plasma radially expand- 
ing in the form of a rarefaction wave [l6| . This plasma 
will be constituted of thermally distributed heavy ions 
with a temperature of the order of the keV (gold) plus 
a component of faster lighter ions from surface contami- 
nants, with an average energy per nucleon of tens of keV 
[l4j,[l5j. Hydrodynamic simulations (carried out with the 
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FIG. 1: Sketch of the typical time of flight experimental ar- 
rangement of the proton imaging technique. Here L = 38 mm 
and I = 4 mm, giving an intrinsic geometrical magnification 
of M w (I + L)/l w 10.5. The long pulse is incident at 45° to 
the target's normal. 



HYADES code [T7|]) indicate that this dense plasma ex- 
pands into a more rarefied background plasma induced by 
the full ionization of the nitrog en g as by secondary X-ray 
emission from the solid target [18( , with an average elec- 
tron temperature and density of 650 eV and 10 16 cm~ 3 , 
respectively. Such an experimental configuration resem- 
bles the one already adopted in Ref . [7| , although with a 
much higher ambient plasma density. The plasma evolu- 
tion was monitored by temporally and spatially resolved 
proton radiographs with a resolution of approximately 
a picosecond and a few microns, respectively [lOj. The 
proton beam was externally generated during the inter- 
action of a secondary laser beam (intensity Il2 ~ 10 19 
Wcm -2 , Short pulse in Fig. [I]) with a 20 /im thick gold 
foil. After having traversed the plasma, the proton beam 
was recorded on a stack of RadioChromic Films (RCF) 
fl9| . The adjustable delay between the two laser beams, 
allowed us to temporally scan the plasma evolution, un- 
til the temporal window in which a collisionless shock 
started to form was found. 

As an example, Fig. [5]a depicts a typical proton ra- 
diograph of the interaction, in which the probing pro- 
ton beam were timed so to traverse the plasma region 
under interest 170 ps after the arrival of the long pulse 
on the gold target. The image shows two main regions 
of proton deflections. A turbulent pattern is visible in 
proximity to the gold target surface (labeled by 'Ablated 
plasma' in Fig. [2Ja) and it is induced by the electrostatic 
and magnetic fields associated with the bulk region of 
the expanding gold plasma (similarly to what discussed 
in Ref. [Tj]). Ahead of this, at ~ 1 mm from the target, 
an approximately circular pattern with a radius of cur- 
vature of ~ 0.9 mm is also visible. As a rule of thumb, 
it is worth recalling that darker grey colors indicate a 
higher proton deposition (proton accumulation) whereas 
lighter grey colors indicate a region of proton depletion. 
Bearing this in mind, this structure is made of two dark 
lines of proton accumulation separated by a lighter stripe 
of proton depletion. Due to the multi-frame capability 
of the proton imaging technique [10| , it has been possi- 




FIG. 2: Two-dimensional proton radiographs: (a) Proton ra- 
diograph of the interaction of a nanosecond laser pulse (red 
arrow) with a 50 /im gold foil (left purple rectangle) corre- 
sponding to 170ps after the arrival of the long pulse (see Fig. 
[l]). Zooms of the proton radiographs of the region highlighted 
by the dashed white rectangle in the frame (a) for different 
proton energies of approximately 11.5 (frame (b)), 10 (frame 
(c)), and 9 (frame (d)) MeV. These energies correspond to a 
probing time of 150, 160 and 170 picoseconds after the start 
of the interaction, respectively. 



ble to follow the temporal evolution of this structure, in 
a single laser shot, within a temporal window of a few 
tens of ps and with a temporal resolution of the order 
of a ps. As an example, Figs. [21b, [5Jc and [5]d depict 
proton radiographs of the same region (highlighted by 
a dashed white rectangle in Fig. [2]a) at 150, 160 and 
170 ps after the arrival of the long pulse, respectively. 
The structure is seen to evolve, while roughly preserving 
its overall shape. Meanwhile, the propagation velocity 
is seen to decrease from an initial value of the order of 
v\ w 2 x 10 6 m/s down to a velocity of approximately 
V2 ~ 8 x 10 5 m/s within a few tens of ps. The ambi- 
ent plasma parameters inferred from HYADES simula- 
tions indicate an ion-acoustic velocity of c s ~ 2 x 10 5 
m/s implying that these velocities would correspond to 
a Mach number of Mi w 10 and Mi ~ 4, respectively. 
We infer that this shock-like structure is collisionless and 
predominantly electrostatic. A collisional shock would 
have a thickness comparable to the ion-ion mean free 
path (i.e. 3 cm for our experimental parameters). A 
magnetized shock would have a thickness of the order of 
the ion-gyroradius. By assuming the nitrogen to be fully 
ionized (as indicated by hydrodynamic simulations), the 
ion-gyroradius can be expressed as: tb ~ rriiV2/(qiB) 
being qi = 7e and rrii = lAm p the nitrogen ion charge 
and mass, respectively. A ion-gyroradius of the order of 
100 microns (measured width of the observed structure) 
would imply an instability-driven magnetic field of the 
order ofB~2 MG, which is unrealistically large in this 
plasma. The electrostatic nature of the shock, which is 
also in agreement with previous experimental [7| and nu- 
merical [16l . l20j work, allows us to quantitatively extract 
the temporal evolution of the electric field distribution 
associated with the shock using Eq. (5) in Ref. [10| (see 
Figs. [2Ia,c)). The electric field evolves from a bell- like 




FIG. 3: The electric field and potential: Experimental electric 
field distribution (a,c) and associated potential (b,d) during 
the transition between the CFDL and the proto-shock. The 
experimental error on the electric field and position are of 
the order of 8% and 5 /im respectively. The profiles refer to 
sequential stages of the same event, and are shown in separate 
graphs for ease of clarity. 



shape (Fig. [3]a) to a symmetric bipolar structure, pass- 
ing through an asymmetrically distributed intermediate 
state (Fig. [He), within 20ps, while keeping an amplitude 
of the order of 80 - 100 MV/m. 

This transition is a clear signature of the precursory 
stages of the formation of a shock, according to the fol- 
lowing dynamics. The propagation of a dense plasma 
into a more dilute one induces a density gradient at 
the blast shell front. Since the thermal speed of elec- 
trons is larger than that of equally hot ions by a factor 
(rrii/m e ) > 40, the electrons diffuse out into the am- 
bient plasma and an ambipolar electrostatic potential is 
set. This structure is called a CFDL [O, [l2J] if the elec- 
tric field is stationary in its rest frame. The progressive 
pile-up of ambient ions by this potential introduces a lo- 
cal density maximum at the blast shell front. Initially 
the density distribution and, thus, the potential on both 
sides of this maximum are asymmetric; the density jump 
is larger towards the ambient plasma. The continuing ion 
pile-up forces the electric field pulse to evolve into a sym- 
metric profile. As soon as the potential is strong enough 
to reflect the incoming ambient and blast shell ions a 
forward and a reverse shock form. We focus on the po- 
tential jump towards the ambient medium. A potential 
jump of approximately 4 kV at £=150 ps [Fig. [2Kb)] can 
not decelerate the ambient nitrogen ions, which move at 
10 c s relative to the blast shell front, even if they are 
fully ionized. The blast shell front is expanding freely. 
The potential [Fig. [3][d)] at t=170 ps has increased to 
12 kV and the structure's speed relative to the ambient 
plasma has decreased to ~ 4c s . This potential is now 
sufficient to reflect ambient nitrogen ions in the reference 
frame of the expanding structure provided their charge 



state is > 4+. The measured Mach number M% ~ 4 
implies that a supercritical forward shock is forming. 

In order to better elucidate the dynamics discussed 
above, a PIC simulation has been performed with ini- 
tial conditions based on those of the ambient medium 
in the experiment. The convergence of the experimen- 
tally measured electric field distribution to a symmetric 
bipolar pulse motivates the study of the collision of two 
symmetric plasma clouds, each with an electron density 
of no = 3.5 x 10 16 cm~ 3 countcrpropagating with a rela- 
tive speed v c ss 5.5c s . The electron and ion temperature 
is set to T — 650 eV. Electrons and ions in each cloud 
are resolved by 4464 computational particles. The left 
plasma cloud occupies the interval — L/2 < x < and 
the right plasma cloud occupies the interval < x < L/2. 
The box length L = 3.6 mm is resolved by 3000 cells and 
we adopted periodic boundary conditions. The plasma 
evolution is followed for T s [ m = 140 ps. 

The intcrpenetration of the two plasma clouds results 
in a density pile-up which sets a net positive electrostatic 
field at the interface. In such an overlap layer, the ions 
are thus slowed down, increasing their local density be- 
yond 2n . Figure SJd reveals at this time a bipolar elec- 
tric field. Its extent (« 100/im) and peak amplitude of 
?s 80 MV/m are practically identical to those observed 
in the experiment (compare Fig. @]d with Fig. [31c). This 
structure develops during 14ps, comparable to the time 
in which the unipolar electric pulse is transformed into 
a bipolar one in the experiment. Later in time (t = 140 
ps), several electric field peaks appear (see Fig. SJe). 
An ion phase space hole [21j is responsible for the sta- 
tionary bipolar electric field at x « 0. Identical shocks 
are located at x ± 60/mi and confine the downstream re- 
gion. We can distrnguish three plasma regions in Fig. 
2)c and x > 0: downstream (x < 60/xm), forcshock 
(60 < x < 600^m) and upstream (x > 600/^m). Two 
unipolar electric pulses, each with a peak electric field of 
~ 0.5 GV/m and propagating at a speed of v s = ±3.4c s 
in the upstream reference frame, are generated. A po- 
tential jump of +5.2 kV at the forward shock as we go 
from the forcshock to the downstream region and a grad- 
ual decrease of the peak forcshock potential from +1.3 
kV at x « 180^tm to the reference value kV at the far 
upstream region are thus set, implying a cumulative po- 
tential between the downstream and upstream region of 
6.5 kV. We have computed these values from an integra- 
tion of the electrostatic field (not shown). This potential 
can reflect the computational ions that move with the 
relative speed 3.4c s and have a mass equal to 250 elec- 
tron masses, generating a beam with density 0.7no (See 
Fig. |4ja). The forcshock potential relative to the up- 
stream accelerates incoming electrons towards the shock. 
The electron core population is at rest at x ~ 1200/im 
and Fig. @]b shows that these electrons have been accel- 
erated to about -5v c at x m 200/im. They gain a kinetic 
energy of about l-2keV. This electron acceleration is en- 
tirely due to the ambipolar electric fields arising from the 
ion density gradient. The potential jump at x ~ 60/im 
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FIG. 4: The simulation results: Simulation plasma at t = 140 
ps. (a) shows the ion phase space density on a linear scale. 
Shocks are located at x±60^m and an ion phase space hole at 
x ~ 0. The electron phase space density is shown on a linear 
scale in (b). (c) shows the distributions of the ion density 
iV(a;) and electron thermal energy density T e (x) normalized 
to those of the upstream, unperturbed plasma, (d) and (e) 
show the electric field profiles close to x ~ at two times. 



enforces a reduction of the typical electron speeds as we 
go from the downstream to the forcshock region, but we 



find fast electrons with v ~ 5v c that leak into the fore- 
shock (Fig. [31b). The electron mixing in the foreshock 
raises the thermal energy density (Fig. 2]c) to about 
1/3 of the downstream one and the ion-acoustic speed is 
raised to c* = \[2c s , where c s is the ion-acoustic speed of 
the unperturbed upstream plasma. The effective Mach 
number of the shock is decreased to M* = 3.4/-\/2 w 2.4. 
This reduction of the effective Mach number might seem 
modest; however, the electron's thermal spread in the 
forcshock depend on the relative speed between the elec- 
trons leaking from the downstream region and the in- 
coming upstream electrons, which have been accelerated 
towards the shock by the potential between the foreshock 
and the far upstream. It is thus independent of the up- 
stream electron temperature. This effect may thus be 
much stronger at SNR shocks, where the downstream 
electron temperature exceeds by far that of the far up- 
stream (ISM) plasma. 

In conclusion, we have reported on the temporally and 
spatially resolved observation of the preliminary stages 
that lead to the generation of a supercritical unmagne- 
tized plasma shock in a laboratory experiment, in condi- 
tions of possible relevance to SNR shocks. A matching 
simulation allowed us to resolve the long-term evolution 
of the structure and it provided insight into the physics 
of the transition layer of this supercritical shock. Both 
experimental data and simulations show that the shock 
potential has the value necessary to reflect some of the 
incoming upstream ions which can in turn efficiently heat 
the foreshock region, thus stabilizing the shock as electro- 
static |22h26| . This effect could reduce the Mach number 
of the forward shocks of SNRs to values well below the 
10 2 — 10 3 that one obtains if one estimates their Mach 
numbers based on the ion acoustic speed of the ISM [27| ■ 
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